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Abstract 

Exact calculations of the radiative corrections to the single spin asymmetry in 
the exclusive pion electroproduction below A(1232) resonance are performed by 
employing the EXCLURAD code. The analysis of the dependence of the obtained 
results on W, 9 cm and <ft shows that absolute values of radiative corrections do not 
exceed 1.5 %. The procedure of sampling the radiative corrections with analytical 
functions has been developed allowing for flexible and rapid applications of the 
corrections in the data analysis. 

1 Introduction 

Evaluation of radiative corrections (RC) to the observables of the exclusive pion electro- 
production 

e" + p — >• e~ + p + 7T° (1) 

rule, performed in the framework of the approach developed in Ref. [1]. Practical 
calculations, which involve application of the code EXCLURAD [2], take a noticeable 
amount of time even in the simple case of several kinematic points. As a result, the 
analysis of the data registered in a wide range of kinematic variables may be hindered 
considerably. Approximations of the results obtained from EXCLURAD with analytical 
functions can substantially decrease time of the data analysis and maintain a high level 
of accuracy in the applied corrections. An input for the approximations is obtained by 
sampling of the energy and angular dependence of RC and the single spin asymmetries 
A and Arc evaluated in the Born approximation and by including all relevant radiative 
processes, respectively. 

2 Kinematics 

The data on single spin asymmetry, which is observed in the reaction (1) in polarized 
electron beam has been collected at MAMI microtron of Mainz [3] . The relevant kinematic 
region is shown in Fig. 1 as a scatter plot with ten points on the W — Q 2 plane. Position 
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Figure 1: Location of the mean values of W and Q 2 corresponding to different experi- 
mental settings in the measurements of the single spin asymmetry at MAMI. 



of the points is determined by the incident beam energy E = 883 MeV, and true average 
of each experimental bin has been found in the analysis of Ref . [4] . 

The code EXCLURAD computes RC to the fourfold cross section d 4 / dW dQ 2 dcos9hd(j)h and 
to polarization beam asymmetry for the reaction (1). The cm system angles 9h and <ph 
are between the direction of the virtual photon momentum (g = k% — k^) and momentum 
of the final pion. As long as a pion in reaction (1) is the undetected hadron, the angles Oh 
and 4>h are defined with respect to the direction of —p', i.e. opposite to the recoil proton 
momentum [1]. In what follows these angles are denoted as cm and (p. 

Calculations of RC, A and Arc have been performed at the fixed value of the virtual 
photon transfered momentum — Q 2 = 0.335 GeV 2 , which corresponded to the central 
value of the Q 2 region represented by the scatter plot in Fig. 1. With the adopted 
definition of <p, both A(<p) and Anc((f>) are symmetric with respect to the angle <ft — 180°. 
In the experiment, cm covered the range 2.5° < cm < 37.5°, and the asymmetries were 
sampled in 10 points inside this interval. The sampling of the <fi dependence has been 
done in 11 points inside in the interval from to 180°. 

As has been demonstrated in Ref. [1], the inclusive asymmetry for the reaction (1) has 
a nontrivial dependence on the center of mass energy W around W ~ 1.3 GeV where it 
changes sign. Therefore, for reliable evaluation of RC one has to sample a much wider 
range of W than that covered by average values in Fig. 1. The relevant values corresponded 
to W = 1.10, 1.15, 1.21, 1.25 and 1.30 GeV. 
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3 Analytical expressions approximating evaluated RC 



3.1 W dependence 

The Bezier curves are used as the basis for approximations of both Born asymmetry 
A ((f), 9, W) and the asymmetry modified by the radiative effects — A RC (0, 9, W). They 
proved to yield a simple approximating tool in the entire range of the kinematics under 
consideration. In this approach the asymmetry evaluated by the EXCLURAD code in a 
number of fixed kinematic points is approximated as 

A^Y.P l B l {u), (2) 

1=0 

where Pi are parameters of the fit and P>i(u) are functions defined on the interval {0, 1}: 

B = (l-u)\ 
B x = 3w(l-w) 2 , 

B 2 = 3u\l-u), (3) 
B 3 = u\ 



Due to this particular choice of functions Bi, Pq and P3 coincide with the values of A on 
the left and right boundary of the interval, respectively Thus, it is only Pi and P2, which 
are to be determined from the fit. 

When fitting W dependence of the asymmetry evaluated on the interval {Wmi n , Wmax}, 
the following change was done: 

W - W min 



u 



w max - w n 



For the applications in the Mainz experiment, u has been calculated by taking W n 
1.10 GeV and W max = 1.30 GeV. 



3.2 Dependence on the azimuthal angle 

Similarly to the case of W dependence, one has to change u = 0/180° in Eq. (2) in the 
fitting of A(4>). On the other hand, one can suggest some simplifications to Eq. (2), which 
would better suit conditions of the Mainz experiment, in which the difference between 
Born and RC asymmetries is very small everywhere except of W ~ 1.3 GeV. This is done 
by noticing that A(<f>) is close to a sine function and by replacing 1st and 4th terms (which 
have to disappear in the case of purely sine function) with a sine term: 

A{4>) = P o sin(0) + P 1 B 1 (u) + P 2 B 2 (u), (4) 

where u = 0/180°. The dependence of parameters Pi on W is obtained by fitting the 
results with Eq. (2). The resulting 2D asymmetry is given by 

A(0,w)=i>f j:Pi m BZ, (5) 

z=o m=0 

where Bq = sin(0), and Bf 2 are usual functions defined in Eq. (3). 
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3.3 Dependence on the angle 9, 



As it is obtained from calculations of the A(9 cm ) by the EXCLURAD, both A(6 cm ) and 
^Rc(^cm) are small and slowly varying functions in the kinematics covered by the exper- 
iment (2.5° < 8 cm < 37.5°). Therefore, the most economical way of approximating the 
asymmetry is the choice of the two-terms functions as follows: 

A{6) = asm(e cm ) + bsm 2 (6 cm ) . (6) 

Correspondingly, approximation of the A(<p, 8 cm ) has been performed for a given value of 
W in the following form: 

A(tf>, e cm ) = £(sin(£ cm )) m ^ PimBi , (7) 

1=0 m=0 

where functions are defined in Eq. (3). Parameters P; m of the fit performed for the 
case with radiative processes are shown in Table 1. 
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m = 1 


m = 2 


m = 3 


1 = 

1 = 1 


-0.51822 
1.3751 


8.2071 
-4.0265 


6.8661 
-6.4145 


-0.51799 
1.3751 



Table 1: Parameters Pi m that determine A RC ((fi,9 cm ) 
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Figure 2: Asymmetry as a function of 6 cm evaluated for three values of the azimuthal angle 
<p — 20° (triangles), 90° (squares) and 160° (circles). Asymmetries calculated in Born 
approximation and by accounting for radiative precesses (assuming AS* = 0.02 GeV 2 ) are 
shown with filled and open symbols, respectively. 
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Figure 3: Asymmetry v4rc(0, #cm) evaluated at W = 1.21 GeV by assuming AS = 0.02 
GeV 2 . 

4 Results 

It was found that the modification of the asymmetry due to radiative effects in the reaction 
under study is the strongest at 9 cm = 37.5°. This is demonstrated by the results of 
calculations A(8 cm ) and A^c(8 cm ) shown in Fig. 2 for three values of the azimuthal angle <fi. 

Radiative corrections depend on the mass window AS* around the missing mass peak used 
in the event selection procedure if AS < ^ max , where z/ max is introduced in EXCLURAD as 
the upper limit of integration in calculating the cross section of the radiative process. The 
parameter z/ max describes the missing mass of the radiative origin (due to undetected 7). 
It is dubbed as "maximum inelasticity" and is defined in Ref. [1] as 

*W ={W- M p f - m\ . (8) 

Undetected gammas additionally distort the missing mass spectrum, which is of the 
hadronic origin (undetected tt° in the considered case). The selection cuts, which are 
usually applied to the distribution of the events versus missing mass squared M%, are 
generally more tight than AS = ^ max , which is modeled in EXCLURAD by replacing the 
upper limit of integration with z/ cut . Accordingly, the mass window AS is related to z/ cut 
as 

AS = z/ max - v cat . (9) 

The corrections naturally increase as the cuts tighten, corresponding to narrowing of the 
mass window AS. 



5 





m = 


m = 1 


m = 2 


m = 3 


1 = 

1 = 1 

1 = 2 


-1.4506 
5.7645 
2.4684 


-4.6186 
12.049 
6.8852 


-1.9493 
7.0872 
8.1888 


0.17369 
-2.7053 
-2.6894 



Table 2: Parameters Pi m that determine A(4>, W) 



The two-dimensional plot for A RC (9 cm ) corresponding to the case of AS = 0.02 GeV 2 
is obtained analytically by employing Eq. (7). It is displayed in Fig. 3. The difference 
between Anc(9 cm ) and A(9 cm ) is better seen from Fig. 4, which shows radiative corrections 
5a defined in Ref. [1] as 

S A = Arc ~ A ■ 100% . (10) 

The set of ten curves in Fig. 4 shows the dependence of 5 a for ten values of 9 cm on the size 
of the AS window used in the event selection procedure. The results have been produced 
for the azimuthal angle = 90°, where corrections are found to be largest. 




0.01 0.02 0.03 0.04 0.05 0.06 



Missing mass window, GeV 2 

Figure 4: Radiative corrections as a function of the AS window evaluated in 10 points 
of the phase space corresponding to the angle 9 cm in the range 2.5° - 37.5° and at the 
azimuthal angle <fi = 90°. 

In this study, a larger number of results have been produced for the kinematics corre- 
sponding to 9 cm = 37.5°. 

The curves in Fig. 5 represent results of the fits, in which Eq. (5) was employed. The 
fitting procedure started from Eq. (4) and three parameters Pi have been obtained in five 
points of W, which allowed to conduct a fit of the W dependence for each of them. It is 
realized by using Eq. (2), and Fig. 6 shows obtained parameters. 

The parameters thus obtained are presented in Tables 2 and 3. They are used in Eq. (5) 
for rapid evaluation of A(<f>,W) and A RC (cf),W) at the angle 9 cm = 37.5° and the AS 
window of 0.020 GeV 2 . The result is displayed in Fig. 7. 
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Fig. 5 Born-term - A, and RC- 
smeared - Arc, asymmetries used 
for the evaluation of radiative cor- 
rections 5a as a function of the az- 
imuthal angle (f) at five values of the 
total cm. energy W - 1.10, 1.15, 
1.21 1.25 and 1.30 GeV. In addition, 
#cm = 37.5° and the AS window is 
0.020 GeV 2 . Full and open points 
correspond to exact calculations of 
the asymmetries A and A RC , respec- 
tively. They are approximated by 
full and dashed lines. 
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m = 1 
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m = 3 


1 = 

1 = 1 
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-1.4544 
5.7668 
2.4694 


-4.6810 
12.365 
7.1451 


-2.2870 
6.2068 
7.6380 


-0.54204 
0.45351 
-0.64024 



Table 3: Parameters Pi m that determine Anc((p,W) 
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Table 4: Parameters Pi m that determine Anc((p, W) 



Dependence of the additive correction defined as A — Arc on and W is shown in Fig. 8 
for the same conditions as in Fig. 7. As it is demonstrated by Fig. 2, the asymmetries 
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Figure 6: W dependence of parameters Pi of Eq. (3) evaluated for Born asymmetry at 
#cm = 37.5°. Three curves (from top to bottom) correspond to Pi, P<i and Pq. 

A(4>, 9 cm ) and Arc(0, 9 cm ) are slow varying functions of 9 cm , which reach maximum at 
the end of the considered kinematic interval - 9 cm = 37.5°. This feature justifies an 
approximate sampling of asymmetries for any given value of 8 cm by scaling down the 
asymmetry found for 9 cm = 37.5°. To this end, one can construct a normalizing function 

}n\4>i @cm) 

A{<j), L) 



>, #cm) 



max 
cm / 
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where A(<f>, 9 cm ) is defined by Eq. (7), and A(4>, 9™ x ) corresponds to the results of the fit 
at 9 cm = 37.5°. The combination of /jv(0, 9 cm ) with the results given by Eq. (5) provides 
a simple and precise tool for evaluating A(4>, W, 9 cm ) in the region around W = 1.2 GeV: 



A(cj>,W,9 cm ) = A( ( j>, W)f N (cf>,9 c 



(12) 
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Figure 7: Single spin Born asymmetry A(<f>, W) (left panel) and the asymmetry Arc(0, W) 
modified by radiative effects (right panel), both evaluated for 9 cm = 37.5°. 




Figure 8: Additive correction obtained as A — Arc, which corresponds to the difference 
between the plots in the left and right panels of Fig. 7. 

Two dimensional plots for A&c{(f),W) evaluated by employing this approximation are 
displayed in Fig. 9. 
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For practical applications the window size AS was varied in RC calculations between 
0.005 and 0.020 GeV 2 . In some plots (e.g. in Fig. 10), for illustration purposes, the 
window size has been increased up to 0.078 GeV 2 . In the kinematics under study (this 
time calculations were done at W = 1.25 GeV), the use of such wide window in the event 
selection serves an educative purpose: it allows to keep virtually all events and obtain 
the asymmetry, which coincides with the Born one. The considered in Fig. 10 case cor- 
responds to rather large phase space of emission of a bremsstrahlung photon. The phase 
space rapidly contracts when W approaches the lower boundary of the kinematic region. 
Therefore, at W = 1.1 GeV all radiatively distorted events can be found inside rather 
small window of 0.005 GeV 2 . On the contrary, as W approaches 1.3 GeV, application of 
the event selection procedure within windows from 0.005 to 0.024 GeV 2 , which is a typical 
size in the data analysis, results in a selection of certain parts of the event samples, which 
have the following properties that can be traced in Fig. 4: 

(1) Large 6 cm — > Large RC — > No dependence on window size 

(2) Small 6 cm — > Small RC — > Certain dependence on window size 

Accordingly, when (AW) 2 approaches its maximum value, the corrections vanish, which is 
seen from converging of curves in Fig. 4 to a single point 5a — 0, and also from convergence 
of measured and Born asymmetries displayed in Fig. 10. 

5 Summary 

In the kinematics of the considered experiment, the radiative effects that manifest them- 
selves in the observed single spin asymmetry are largest at W — 1.3 GeV and at the 9 cm = 
37.5°. But even there they do not change the Born asymmetry by more than 1.0 - 1.5 % 
(in absolute value). A monotonous rise of A{6 cm ) and A&c{6 cm ) in the angular interval (0, 
37.5°) allows to employ simple analytical approximation of the asymmetries and radia- 
tive corrections by evaluating them in the maximum (0™ x = 37.5°) and scaling down to 
any value of 9 cm inside the considered interval. A number of analytical functions, which 
approximate numerical results from EXCLURAD, offer flexible tools for fast applications 
of radiative corrections as a function of W, 9 cm and (f>. 

The dependence of radiative corrections on the AS window, if it is chosen as 0.005 < AS < 
0.024 GeV 2 , can be neglected in the whole range of W but from the very different grounds: 

(1) at W = 1.1 GeV any AS window completely covers the phase space of radiative events, 

(2) at W — 1.3 GeV any AS window is small compared to the phase space of radiative 
events, but falls onto the plateau of the dependence of the radiative corrections on the 
AS window size. 
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Figure 9: Single spin asymmetry ^4 R c(0, W), which includes radiative effects evaluated 
for the values of the polar angle 9 cm = 2.5° (upper left), 13.8° (upper right), 30° (lower 
left) and 37.5° (lower right). 
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Figure 10: Single spin asymmetry as a function of the AS" window evaluated at two data 
points corresponding to the angles 6 cm = 18° and 37.5°. Full and dashed lines show A and 
Arc, which are the Born term and smeared by radiative effects asymmetries, respectively. 



7 Appendix 

This section explains how the Fortran functions written for the approximation of the 
single spin asymmetries and radiative corrections can be used in practical work. 

All asymmetries approximated by the functions correspond to the MAID-2000 parameters 
used as the input for EXCLURAD. 

The arguments of the functions are in units of "GeV" (for the total cm. energy W) and 
"degree" (for angles 9 cm and 0). 

1. FBORNA(PHI,W) 

Returns asymmetry A((p, W) (%) in Born approximation at fixed 8 cm = 37.5°. 
It has been used in producing left-panel plot of Fig. 7. 

2. FRC020(PHI,W) 

Returns asymmetry A RC (0, W) (%), which includes radiative effects, at fixed 9 cm 
= 37.5° and corresponds to the event selection in the missing mass window of 0.02 
GeV 2 . 

It has been used in producing right-panel plot of Fig. 7 as well as lower-right plot 
in Fig. 9. 
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3. VRADCOR(PHI,W) 

Returns correction obtained by subtracting asymmetry A^c(<f), W) from Born asym- 
metry A(<f>,W). Uses functions FBORNA(PHI,W) and FRC020(PHI,W). 
It has been used in producing the plot in Fig. 8. 

4. FASYMT(PHI,TH) 

Returns Anc((f>, 6) (%) at fixed W = 1.21 GeV and corresponds to the event selection 

in the missing mass window of 0.02 GeV 2 . 

It has been used in producing the plot in Fig. 3. 

5. FNORM(PHI,TH) 

Returns the ratio of Arc(</>, #cm) and Ajic(4>, &™ x ) f° r an Y given 9 cm in the W region 
around the mean value of the considered experiment - W « 1.2 GeV. It has the same 
structure and uses the same set of fit parameters as the function FASYMT(PHI,TH). 
It has been used in producing the plots in Fig. 9. 
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